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We address the limit ol the Gaussian operations and classical communication in the problem of 
quantum state discrimination. We show that the optimal Gaussian strategy for the discrimination 
of the binary phase shift keyed (BPSK) coherent signal is a simple homodyne detection. We also 
propose practical near-optimal quantum receivers that beat the BPSK homodyne limit in all areas 
of the signal power. Our scheme is simple and does not require realtime electrical feedback. 

PACS numbers: 03.67.Hk, 42.50.Dv 



I. INTRODUCTION 

Discrimination of the binary phase shift keyed (BPSK) 
coherent states {|q:),|— a)} with the minimum error is 
one of the most fundamental issues in optical commu- 
nication and quantum signal detection theory. Coherent 
communication theory has been developed based on semi- 
classical theory where these signals are detected by ho- 
modyne measurement. For the signals with equal prior 
probabilities, the average error probability is given by 
Perr — erfc[-\/2| a |] /2. This is often called the shot noise 
limit or the homodyne limit, falling short of the conven- 
tional error free criterion (10~^) when jap < 10. 

It is, however, well known that the quantum optimal 
receiver can largely surpass the homodyne limit. The 
optimal measurement is mathematically given by a two- 
dimensional projection measureme nt and it at tains the 
minimum error probabihty of (1 - Vl - e-^l"!") /2 which 
is called the Helstrom bound Ij. Kennedy proposed a 
simple near-optimal receiver using a coherent local os- 
cillator (LO) and photon counting Its error rate is 
only twice larger than the Helstrom bound and is smaller 
than the homodyne limit when |ap > 0.4. Dolinar 
then extended this 'Kennedy receiver' to the optimal one 
by introducing the adaptive electrical feedback which is 
enough faster than the optical signal pulse width 3)] (see 
also [l|, 0, H, H, 01)- Although the Dolinar's concept has 
been demonstrated recently [1], it is still challenging to 
experimentally beat the homodyne limit with this ap- 
proach because of its complicated system. 

From a quantum mechanical point of view, homodyne 
measurement belongs to the class of Gaussian operations, 
i.e. described by up to the second order nonlinearity, 
while photon counting is the non-Gaussian one. The role 
of Gaussian operations in quantum information protocols 
Q or quantum state estimation [10] have been widely in- 
vestigated. On the other hand, it has also been shown 
that some of the important protocols cannot be per- 
formed by only Gaussian operations and classical commu- 
nication (GOCC) and inevitably requires non-Gaussian 
operations, e.g. quantum computing (ll| . entanglement 



distillation of Gaussian states 12 , [l^ [13] , and the opt 



mal cloning of coherent states 

In this paper, we first show that the homodyne mea- 
surement is the best strategy to discriminate the binary 
coherent states within GOCC. To our knowledge, this is 
the first result addressing the Gaussian limit in quantum 
state discrimination scenario. In the second half of the 
paper, we propose novel non-Gaussian quantum receivers 
based on the Kennedy receiver, that beat the Gaussian 
limit for any ja]^. In particular, we point out that the 
amount of displacement in the Kennedy receiver is not 
optimal. Our schemes do not require realtime feedback 
and are simple and practical to experimentally overcome 
the homodyne limit with current technology. 



II. DISCRIMINATION VIA GAUSSIAN 
OPERATIONS AND CONDITIONAL DYNAMICS 

In this section, we show that "the minimum error dis- 
crimination of a set of two coherent states {ja},]— a)} 
with the prior probabilities {p4^,p^} under GOCC is at- 
tained by the homodyne detection" . For simplicity, a is 
assumed to be real. Gaussian operation is defined as the 
operation that maps Gaussian states to Gaussian states. 
For Gaussian input states, properties of Gaussian opera- 
tions have been well investigated [13, [3] • In our problem, 
however, although each signal state is given by a Gaus- 
sian state, the signal from the receiver's viewpoint is an 
ensemble of these states, pi — p+la)(a) -\- p-\—a){—a\, 
i.e. non-Gaussian. This is because the receiver does not 
know which state he or she is receiving. We therefore 
start by revisiting the measurement processes based on 
GOCC. 



A. Characterization of the measurements with 
GOCC 

It is known that any Gaussian operation (Gaussian 
completely positive (CP) map) in optical system can be 
implemented by adding an ancillary system prepared in 
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Gaussian state, applying Gaussian unitary operation on 
the whole system (implementable via linear optics and 
squeezing), and then discarding and/or performing ho- 
modyne measurements on a part of the system fl^ . The 
CP map including measurements is not always trace- 
preserving and the output quantum state might be condi- 
tioned on the measurement outcome. When Gaussian op- 
erations are sequentially applied, the measurement out- 
comes (classical information) are sometimes useful to dy- 
namically renew each step of quantum operations, which 
is called conditional dynamics. In the following, we char- 
acterize two types of measurements consisting of Gaus- 
sian operation with and without conditional dynamics. 

The first one is the measurement with only Gaus- 
sian quantum operation (without conditional dynam- 
ics). Here we call it a 'Gaussian measurement'. A 
generic physical model of the Gaussian measurement is 
depicted in Fig.[lja), consisting of Gaussian unitary op- 
eration, Gaussian ancillary states, and homodyne detec- 
tors. After the Gaussian measurement is finished, a set 
of measurement outcome is classically post-processed, in 
our case, to make the decision which signal is detected. 
Throughout this section, we consider only 'noise-free' op- 
erations and measurements, that is we assume that an- 
cillary states are always pure and the system is never dis- 
carded. It does not lose generality. Since to discard some 
of the modes means to lose their information, it is realized 
by measuring them via homodyne detectors and ignor- 
ing the outcomes, where the latter is included in classical 
post-processing. Mixed ancillary states are provided by 
first preparing entangled pure states and then discarding 
a part of them, therefore, can be prepared by pure states 
and post-processing. These imply the generality of the 
noise-free model. A Gaussian measurement detecting an 
L-mode quantum state is mathematically described by a 
positive operator-valued measure (POVM) {nG(r, (5)}5 
where nG(r,(5) > 0, / d6tlG{T,S) = /, and / is an 
identity operator (see Appendix for its derivation from 
the physical model). The operator IlG(r,(5) is a Gaus- 
sian operator, i.e. its characteristic function is described 
by x(a;) = exp[— iw-^Fw -I- iS'^uj] where F and 5 are the 
2L X 2L covariance matrix and the 2L-dimensional dis- 
placement vector, respectively, and T is the transpose 
operation. 

The second one consists of GOCC which we call a 
'GOCC-measurement'. As illustrated in Fig. [T{b), its 
generic model is described by sequential Gaussian oper- 
ations with conditional dynamics via classical communi- 
cation. Each step of Gaussian operation includes Gaus- 
sian operation and ancillary states, and a Gaussian mea- 
surement detecting a part of the system. The measure- 
ment outcomes are applied to modify the following step 
of Gaussian operations in realtime, which is the condi- 
tional dynamics via classical communication. After the 
whole quantum measurement process is finished, all of 
the measurement outcomes is used for the classical post- 
processing. The whole process except the post-processing 
is described by a POVM {nGOCc(F(a;), (5(a;))}^ with a 
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FIG. 1: (Color online) Generic physical models of (a) Gaus- 
sian measurement and (b) GOCC-measurement. The solid 
and dotted lines represent quantum and classical signals, re- 
spectively. HD: homodyne detector, GU: Gaussian unitary 
operation, GM: Gaussian measurement, PP: post-processing. 



covariance matrix T{x) and a displacement S(x) where 
nGOCc(F(a;),<5(a;)) > and / dxIlGOCciTix),6{x)) ^ 1. 
The parameter x indicates the pattern of which con- 
ditional dynamics is applied during the whole process. 
Again, throughout the section, we restrict the GOGG- 
measurements to be noise-free. 

Gaussian measurement is well characterized and eas- 
ily applied to the optimization problem on the state 
discrimination. On the other hand, although GOCC- 
measurement is also well defined, it is not easy to handle 
its POVM directly. To prove the optimality of the ho- 
modyne limit under GOGC-measurements, therefore, we 
first show that the optimal Gaussian measurement with- 
out CC is a homodyne measurement. Then we discuss 
an important property of the conditional output from 
a Gaussian operation with an input of binary coherent 
state signals. Finally, we prove that even in the GOCC- 
measurement scenario, conditional dynamics is not useful 
and thus a simple homodyne measurement is optimal. 



B. Optimal Gaussian measurement 

A Gaussian measurement for the single-mode input 
state is described by {II{Tm, Dm)}dm where Dm is a 
two-dimensional vector and 



c_ s 
s c+ 



(1) 



where c± = cosh(2r) ± sinh(2r') cos ip, s = sinh(2r) sin (p, 
and r and if are the real parameters. The mini- 
mum error probability to discriminate {|±Q;),p±} by 

{II(Tm, Dm)}dm with given r and (p can be calculated 
from the probability distribution to detect each signal 
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P±{Dm) = (±a|n(r7K, D7n)|±a). Applying a conven- 
tional Bayesian decision strategy as a post-processing, we 
obtain 



p(G) 
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where 



1 + cosh(2r) + sinh(2r) cos cp 
2(1 -Fcosli(2r)) ■ 



(2) 



(3) 



(G) 

It is apparent that Pe is minimum when ip = and 
r = oo, which implies that the homodyne detection with 
the phase iy9 = is the optimal strategy within all possible 
Gaussian measurements. 



C. Conditional output states from a Gaussian 
operation 

As mentioned above, a GOCC-measurement consists 
of a sequence of Gaussian operations that include partial 
measurements. In this subsection, before discussing a 
whole GOCC-measurement process, we pick up one step 
of the sequence and address a useful property of the con- 
ditional output from a Gaussian operation with the bi- 
nary coherent state inputs. 

Let us consider the noise-free conditional Gaussian 
operation which transforms a single-mode input to an 
A'^-modc output where the output state is conditioned 
on a partial measurement outcome (Im- Suppose an 
input state is |a) or | — a). We show that for any 
d^K, the conditional output states of the inputs | ± a) 
can always be transformed to |Q;'j_)(a'^| ® p'^^^ by the 
same -independent deterministic Gaussian operation, 
where a'j. = zta' -I- a'{dM) f^nd a' is independent of 
cIm- We also show that, for a statistical ensemble input 
Pi — p+\a){a\ + p-\ — a){~a\, the conditional output is 

where 



similarly transformed to be p'^ p^^^ 

p',^p'+{dM)\a'+){a'^\ +p'_{dM)\a'_){a'_ 



(4) 



Let /5(7, d) be a density matrix of a Gaussian state with 
a covariance matrix 7 and a displacement d. For exam- 
ple, the coherent states |±q;) are denoted as /5(/2,±da) 
where ±da = [±\/2a, 0]^ and I2L is a 2L x 2L iden- 
tity matrix. The conditional operation is described as 
follows. The initial single-mode state is interacted with 
M — 1 Gaussian auxiliary states (M > N) via a Gaussian 
unitary operation. Without loss of generality, we can set 
the auxiliary states to be M — 1 vacua. At the covariance 
matrix level, Gaussian unitary operation is described by 
the matrix transformation via a symplectic matrix S and 
an additional displacement d. These transform the coher- 



ent states /5(/2,±da) as 
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(5) 
(6) 



where S and 7 are 2M x 2M matrices and d and d are 
2Af -dimensional vectors. The A^-mode conditional out- 
put is obtained by performing an (Af — A^)-mode noise- 
free Gaussian measurement {11(7^, dx)}^^ . For con- 
venience, we divide the system by the first N modes and 
the remaining M — N modes and call them the system 
A and B, respectively, as 



(7) 



After performing the Gaussian measurement on the sys- 
tem B, each of the signals (| ± a)) is transformed to an 
A'^-mode conditional output state as 
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Note that each of p{Tout, D±) is a pure state since the 
operations are noise-free. 

Let us show that /5(rout,£'±) can be simultaneously 
transformed to |a'^)(aj_| (8)paux via Gaussian unitary op- 
erations. Since each of p{Tont, D±) is a pure state, there 
exists a symplectic transformation (i.e. Gaussian unitary 
operation) So such that [l3| 



Tout SD^outSjj — I2N, (11) 

where the displacement is also transformed as 



±SdD + SdDm- 



(12) 



Note that Sd depends only on Tout and thus inde- 
pendent of dM- Let SdD = [di,d2,--- ,d2N]'^ and 
SdDm = [di,d2,--- ,d2N]'^- We can transform them 
to [±d', 0, • • • , 0]^ and [d[,d2,- ■ ■ ,d2j^]'^ , respectively, 
by some combination of linear optics (beamsplitters and 
phase shifters) where the covariance matrix I2N is kept 
to be invariant. Again parameters of the beamsplitters 
depend only on {di}i, and independent of {di}i, i.e. free 
from dM- After these operations, the states are trans- 
formed to be the desired ones 



p{^out,D±) p{l2, [±d' + d{,d'^f 



(13) 
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where p'^^^ is a product of iV — 1 coherent states with the 
displacement [d'^.d'^,- ■ ■ .d^^^]^- These are the desired 
ones. 

Finally we apply the above scenario onto the initial 
state of Pi — p+p{l2,+da) + p-p{l2,—da)- Following 
the above procedures, its conditional output after the 
Gaussian operation is given by 



Pout ==p+P+(d7n)/5(rout,-D+) + p- P-{dM) p(r out -.D^), 



where 
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(14) 



(15) 



After the unitary operation of Sd and appropriate linear 
operations, the state is transformed to be 



Pout ^Pi®Pe. 



where 



p+P+{dM)p{hAd' + d\,d'2Y) 
+P-P^{dM)p{l2A-d' + d'i,d'2_ 



D. Optimal GOCC- measurement 
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FIG. 2: (Color online) Schematic of the near-optimal quan- 
tum receivers, (a) Type-L Photon detector -I- optimal Gaus- 
sian unitary operation, (b) Type-II: Photon detector -|- opti- 
mal displacement. 



transformation pout ~^ Pi and the homodyne detection 
where any parameters in those processes are indepen- 
dent of djVi. It implies that the conditional dynamics 
is not necessary for to optimize the second step mea- 
surement [l8l |. Consequently, the optimal whole process 
of these steps is described by a Gaussian measurement 
and thus, as already shown, is a homodyne measure- 
ment. Note that this statement is obtainable without 
specifying a concrete process of the first step Gaussian 
operation. An extension of the above scenario to the 
multi-step one is straightforward, which proves the opti- 
mality of the homodyne measurement within all possible 
GOCC-measurement. 



III. NEAR-OPTIMAL QUANTUM RECEIVER 
BY USING A PHOTON COUNTER 



Let us now turn to the state discrimination via GOCC- 
measurements. To specify the role of conditional dy- 
namics, we first consider a simpler measurement sce- 
nario where a single conditional Gaussian operation and 
a Gaussian measurement are sequentially operated on the 
signal. Denote the partial measurement outcome at the 
former step as dj^, which is informed to the latter mea- 
surement step to optimize the process of Gaussian mea- 
surement. After these GOCC processes, all measurement 
outcomes are classically post-processed. 

Applying the result in the previous subsection to the 
initial state of pi, the conditional output from the first 
Gaussian operation can be transformed to correspond- 
ing p'^ described in Eq. p?]) via -independent de- 
terministic Gaussian operations. Let this operation be 
a part of the second step Gaussian measurement (if 
necesssary, one can add pf^y^ as an ancilla). Then the 
remaining task in the measurement is to discriminate 
two coherent states {|a+),|a'_)} with the prior prob- 
abilities of {p'-^{dM),p'-idM)}- As already mentioned, 
the optimal Gaussian measurement is given by a sim- 
ple homodyne detection. Its phase (p is determined by 
geometric configuration between a'_^_ and a'_ and since 
a± = {±d' + d[+ i4)/v^ (see Eq. it is always 

given by (y3 = which is irrespective to the values of d'^ , 
d'2 and thus also dM- An optimal strategy for the sec- 
ond step Gaussian measurement therefore consists of the 



The homodyne limit (GOCC limit) stated in the pre- 
vious section is overcome by adding a non-Gaussian 
measurement device. In this section, by extending the 
Kennedy receiver, we propose a simple near-optimal re- 
ceiver where a photon counter, which is a typical non- 
Gaussian operation device, is added to the prior Gaussian 
operation. In what follows, we assume = p- = 1/2 
for simplicity. 

In the Kennedy receiver, the BPSK signal {|a), | — a)} 
is shifted to {|2a),|0)} by the displacement operation 
D{a) — exp(aat — a*d), where a and d^ are annihilation 
and creation operators, respectively, and then detected 
by an on/off type photon detector which discriminates 
zero or non-zero photons. It is well known that D{a) 
can be realized by using a beamsplitter with the trans- 
mittance r ^ 1 and the coherent LO \a/^/T). Here, 
we extend the Kennedy receiver and consider the setup 
depicted in Fig. HUa), where the displacement D{a) is 
replaced by a single-mode general Gaussian unitary op- 
eration Ug- We will seek its optimal structure. 

An on/off detector is described by the POVM 
{noff,non} with 

00 

floff = e-" ^ (1 - vr\m){m\, Hon = / - ftofi, (18) 
where |m) is an m-photon state, rj is the quantum effi- 
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ciency, and v is the dark counts. The Gaussian unitary 
operation JJq consists of phase shift, displacement, and 
squeezing, while one can omit the phase shift since the 
on/off detection is insensitive to the global phase. Then 
the average error probability is calculated from 



Pe = 



[a\U:^U,sUG\a) 



'(19) 

where Ua = DiP)S{C), 5(C) = exp[i(C*a2-Ca^2)] is the 
squeezing operator, and C. — re'f is the complex squeez- 
ing parameter. 

After some algebra, one can find that r, ip, and (3 have 
the extreme at the same point, where (p = 0, f3 is real, 
and the optimal displacement /3opt and squeezing Topt, 
are given by /3 and r satisfying 



SrjaP _ ( A-q{a'^ + (3"^) ?7 + (2 - T])e' 



□4r 



\ — -e 
X tanh 



2 '-^2) + (2 - 7y)e-2' - 
4'- ?? + (2 - ?7)e2'- 



a = (3 tanh 



7/ + (2 - r;)e-2'- 

Arjal3 
77 + (2 — ■q)e 



-2r 



(20) 
(21) 



simultaneously. The optimized average error probability 
is then given by 



+ (2 - 77)e2'-°pt)(?7 + (2 - 77)e-2'-opt) 



X exp 
X sinh 



2y?(a2 + /3gpJ 
77 + (2 - 77)e-2''°pt 

4?7a/3opt 



r/ + (2 — ?7)e ^'"opt 



(22) 



In the following, we call it as the Type-I receiver. It 
should be noted that if one can use an arbitrarily higher 
order nonlinear unitary operation instead of J7gj it is able 
to achieve the Helstrom bound rigorously ^ . 

On the other hand, if one is restricted to use only linear 
unitary operation, that is the displacement -D(7), the 
conditions in Eqs. (pOI and (PT|) are simplified as 



7 tanh(2?7Q!7), 



(23) 



The schematic is shown in Fig. [^fb) and we call it the 
Type-II receiver. Its average error probability is given by 



2 



p'^ sinh(2r/a7opt) , 



(24) 



where 7opt is the 7 satisfying Eq. p3|) . Note that its 
physical setup is the same as that of the Kennedy re- 
ceiver. However, we stress that 7opt 7^ a in general 
and thus the conventional Kennedy receiver is easily im- 
proved by using -D(7opt) instead of D{a). Figure [3l^a) 
plots the average error probabilities for the Type-I, Type- 
II, and Kennedy receivers, the homodyne limit, and the 
Helstrom bound while Topt, /3opt, and 7opt are shown in 
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FIG. 3: (Color online) (a) Average error probabilities for the 
Type-I, Type-II, and Kennedy receivers, the homodyne limit, 
and the Helstrom bound, (b) The optimal displacements and 
squeezing for the Type-I and II receivers. 



Fig. |3)^b) . It is shown that the error probabilities for both 
the Type-I and II receivers are better than the homodyne 
limit for any |a|2. 

Let us finally discuss the practical perspective of 
these non-Gaussian receivers, particularly, the Type-II 
receiver. The superiority of the Type-II receiver rather 
than the Kennedy receiver in jap < 1 is significant to 
beat the homodyne limit in realistic experiments. It 
is known that the Kennedy-type receiver is not robust 
against thermal noise or dark counts [l^. Moreover, even 
without environmental noises, the mode mismatch be- 
tween the signal and LO causes additional dark counts. 
As mentioned above, the displacement D{P) is realized 
by interfering the signal with the coherent state LO 
l/3/vl — t) via the beamsplitter of the transmittance r. 
The effect of mode mismatch can be characterized by 
introducing the mode match factor ^ (0 < ^ < 1) rep- 
resenting the overlap between the signal and LO pulse 
areas. Since these two pulses are in a coherent state, the 
average intensity of the signal field after the interference 
is simply given by 



/=(l-e) (T|a|2 + |/3|2)+^|±V7a + /3| 



(25) 



Due to its Poissonian photon number distribution, the 
average discrimination error including r and f at the 
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FIG. 4: (Color online) Average error probabilities for the ideal 
homodyne receiver and the Type-II and Kennedy receivers 
with practical imperfections, r — 0.99, rj — 0.9, u — 10~^, 
and C = 0.995. 



the homodyne limit for any signal photon number re- 
gion. Because of the recent experimental progress of high 
efficiency photon detectors [2l| and universal squeezing 
operations j2^ . we believe that now it would be feasi- 
ble in near future to beat the homodyne limit in digital 
optical communication experiments. 

Note added: Proof-of-principle experiment of the 
Type-II receiver is recently demonstrated [23i] . 
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on/off detector is described as 

pF^I~ e-''"''(™'+^-') sinh (277^V^a7opt) , (26) 
where 7opt fulfills the optimality condition 

^^/ra = 7opt tanh(2ry^a7opt)- (27) 

An example of the average error probabilities including 
the imperfections is shown in Fig. 0] which clearly shows 
the advantage of our proposed receiver would be crucial 
to experimentally observe the gain of the non-Gaussian 
measurement beyond the homodyne limit. Although the 
requirement for rj in the weaker signal is still high, recent 
experimental progress in this field is rather promising 



IV. CONCLUSIONS 

In this paper, we have addressed the discrimination of 
the BPSK signals and proved that the homodyne limit 
is the minimum error probability attainable via Gaus- 
sian operations and classical communication. This is the 
first clarification of the limit of Gaussian operation in the 
state discrimination problem. Although it is shown for 
the binary coherent states that any conditional dynam- 
ics due to CC is not effective, we note that this would 
not be the case for the discrimination of more than two 
signals. Related to this topic, an increase of the mu- 
tual information by the adaptive homodyne strategies 
has been numerically observed [20|. For further investi- 
gation into this direction, more useful formulation of the 
GOCC-measurement would be necessary. 

We have also proposed the near-optimal quantum re- 
ceivers for the BPSK coherent signals, that are based on a 
photon detector and Gaussian operations. Our schemes 
are simple and do not require realtime electrical feed- 
back although their error probabilities are better than 



APPENDIX: GENERIC MODEL OF THE GPOVM 

Here we show that the POVM of the physical model 
illustrated in Fig.[TJa) is always described by {HgCF, d)}s 
introduced in Sec. IIA. Let and pf^^ be an iV^i-mode 
input state and an A^^-mode ancillary state, respectively. 
The probability distribution of the measurement is given 
by 



p(c?hd) 



Tr^B 
TrA 



AB 



(A.l) 



where Ug^ is an {Na + NB)-uiode Gaussian unitary op- 
eration and {n^Q(rHD, c^hd)} represents Na + Nb ho- 
modyne detectors with the measurement outcomes de- 
noted by duD • Note that homodyne detection is a Gaus- 
sian measurement (projection onto an infinitely squeezed 
states) and thus characterized by the covariance matrix. 
In Eq. dMl), Fhd is a 2(7Va + A^b) x 2(iVA-)-7VB) diagonal 
matrix and 



HD 



diag[e 



-2r 2r ~2r 2r 



(A.2) 



with r — > oo. 

Equation (|A.1|) imphcs that the POVM of the Gaus- 
sian measurement model is derived from a set of op- 

erators {TrB[/3f,,J7^^tn^g(rHD, dHD)C>#^]}dHD- Let 
us describe it by covariance matrices and displace- 
ments. Denoting the symplectic transformation cor- 
responding to Ug^ by 5, the unitary transformation 



HD,rfHD)C^S^ described by 



Fhd 

C^HD 



Fs, 



Sd 



HD 



is- 



(A.3) 
(A.4) 
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Then after tracing out the ancillary system B, we find 
that the above POVM is given by an A^^^-mode Gaussian 
operator IIg{T,6) with 



r = r^-rc 

S = dA-Tc 



raux + 
1 



B 



(A.5) 
(A.6) 



where we have denoted 
Ts = 
ds = 



Ta Tc 

dA 
ds 



(A.7) 
(A.8) 
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